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ABSTRACT 

Surveys of distant galaxies with the Hubble Space Telescope and from the ground have shown that 
there is only mild evolution in the relationship between radial size and stellar mass for galactic disks 
from z ~ 1 to the present day. Using a sample of nearby disk-dominated galaxies from the Sloan Digital 
Sky Survey (SDSS), and high redshift data from the GEMS (Galaxy Evolution from Morphology and 
SEDs) survey, we investigate whether this result is consistent with theoretical expectations within the 
hierarchical paradigm of structure formation. The relationship between virial radius and mass for dark 
matter halos in the ACDM model evolves by about a factor of two over this interval. However, N-body 
simulations have shown that halos of a given mass have less centrally concentrated mass profiles at 
high redshift. When we compute the expected disk size-stellar mass distribution, accounting for this 
evolution in the internal structure of dark matter halos and the adiabatic contraction of the dark 
matter by the self-gravity of the collapsing baryons, we find that the predicted evolution in the mean 
size at fixed stellar mass since z ~ 1 is about 15-20 percent, in good agreement with the observational 
constraints from GEMS. At redshift z ~ 2, the model predicts that disks at fixed stellar mass were 
on average only 60% as large as they are today. Similarly, we predict that the rotation velocity at 
a given stellar mass (essentially the zero-point of the TuUy-Fisher relation) is only about 10 percent 
larger at z ~ 1 (20 percent at z '--^ 2) than at the present day. 

Subject headings: galaxies: spiral - galaxies: evolution - galaxies: high redshift - surveys - cosmology: 
observations 



1. INTRODUCTION 

The relationship between the radial size and the lu- 
minosity or stellar mass of galactic disks is a funda- 
mental scaling relation that reveals important aspects 
of the formation history of these objects. The size, lu- 
minosity/mass, and rotation velocity for m a 'fundamen- 
tal p lane' for disks at the present epoch (jPizagno et al.l 
l2006f) that is analogous to the more familiar fundamental 
plane for early- type galaxies (Burstein et al. 1997). The 
zero-point, slope, and scatter of the fundamental plane 
for both disks and spheroids, and the evolution of these 
quantities over cosmic time, pose strong constraints on 
models of galaxy formation. 



In the modern paradigm of disk formation. Cold Dark 
Matter (CDM) dominates the initial gravitational po- 
tential, and dark matter (DM) and gas aquire angu- 
lar mornenturn via tidal torques in the early universe 
(jPeeblesI Il969t) . When the gas cools and condenses, 
this angular momentum may eventually halt the col- 
lapse and l ead to the formation o f a rotationally sup- 
ported disk ([^nr&^^fetathio^IiH^)- Under the assump- 
tion that the specific angular momentum of the pre- 
collapse gas is similar to that of the DM and is mostly 
conserved during collapse, this picture leads to predic- 
tions of present-day disk sizes that are in reasonably 
good agreement with observations (iKau ffmannI '199^ 
Dalcanton et al.lll997l: iMo et al.lfl99l:lAvi la- Reese et a! 
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19981; ISomerville & Primack"1999'; 'v an den Bosch i2000l : 
Dutton ct al. 2007; Gncdin ct al. 200^^ 

However, in the most sophisticated numerical hydro- 
dynamic simulations of disk formation in a CDM uni- 
verse, the proto-disk gas tends to lose a large frac- 
tion of its angular momentum via mergers, leading to 
disks that are too small and comp act (Navarro & Whit3 
'1994; Som mer-Larsen et al.n i999: N avarro fc Steinmct^ 
|200Q,). It is still not clear whether this problem re- 
flects a fundamental difficulty with CDM (i.e., excess 
small scale power), or is due to inadequate numeri- 
cal resolution or treatment of "ga strophysical" processes 
like s tar formation an d feedback l|Governato et al.ll2004l : 
R obertson et al.l [20041 ). However, it has been suggested 
that delayed cooling and star formation, perhaps due 
to strong feedback in low-mass proge nitors, could help 
to stem this angular momentum loss (jWeil et al.l 119981 : 
iMaller fc Deke]|l2002f ). These ideas can be tested by ob- 
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serving the evolution of disk scaling relations at high red- 
shift relative to the present epoch. 

The observational relationship between radial size (ef- 
fective radius or disk scale length) and luminosity or 
stellar mass for disks at low redshift has now been well- 
characterized by s tudies based on t he Sloan Digital Sky 
Survey (SDSS; e.g. lShen et alEoO^ hereafter S03). Sev- 
eral pioneering studies in the past decade studied the 
size-l uminosity relation for disks out to r edshift z ~ 1 
(e.g. iLillv et all flOOSt iSimard et all [l999l) . Sizes were 
also measured f or Lyman Break Galaxies (LBGs) at 
redsh ifts z ~ 3 (jGiavalisco et al.l 119961 : iLowenthal et al.l 
Il997( ). However, there was significant disagreement be- 
tween the results of different studies, and a clear under- 
standing of the evolution of galaxy sizes over cosmic time 
was hampered by the difficulty of obtaining large samples 
with accurate redshifts and high-resolution imaging, and 
by concerns about the impact of surface brightness se- 
lection effects. An additional problem with interpreting 
high redshift data is that as higher redshifts are probed, 
the observed optical begins to shift into the rest-UV, and 
k-corrections become highly uncertain. Moreover, the 
stellar mass-to-light ratios of galaxies increase as their 
stellar populations age, and therefore the size-luminosity 
relation would change over time even for galaxies that 
were just "passively" aging. 

Recently, new studies with the Advanced Cam- 
era for Surveys (ACS) on the Hubble Space Tele- 
scope (HST) have provided greatly improved constraints 
on the disk size-lu r nino sity relati on at high r e dshift . 
iRavindranath et all (|20b 4) and F erguson et all (|2004D 
presented size distributions for 2: < 1 disk-type galaxies 
and for rest-UV selected galaxies from 1.4 ^ z < 6, re- 
spectively, based on samples selected from the Great Ob- 
servatories Origins Deep Survey fGOODSl.'B arden et all 
(jlOOS, hereafter BOS) presented the luminosity-size and 
stellar mass-size relation out to 2; '--^ 1 based on the 
GEMS (Galaxy Evolut ion from Morphology and SEDs) 
survey (|Rix et al.ll2004D . They concluded that disk galax- 
ies of a given size at 2: ~ 1 are ~1 magnitude brighter in 
the F-band, but that there is less than about a ten per- 
cent change in the stellar mass at a given size between 
z ~ 1 and the present. This result is consistent with a 
mean stellar mass-to-light ratio that increases with time, 
as expected based on the simple aging of stellar popula- 
tions. An interesting co mplementary result was recently 
found by iSargent et all (j2006|) based on the COSMOS 
survey. They found that the number density of disks 
with half-light radii between 5 and 7 kpc is nearly con- 
stant from z ~ 1 until the present. Coupled with findings 
that the stellar mass function of spir al galaxies at z ^ 1 
was about the same as it is today (jBorch et al.ll2006l ). 
this reinforces a picture in which disks have a fixed rela- 
tionship between their stellar mass and their radial size 
over this redshift interval (i.e., that as disks grow in mass 
through star formation, their sizes grow in such a way as 
to kee p them on this re lation, on average). 

iTruiillo et al.l (|2004f ) measured luminosity-size and 
stellar mass-size relations in the rest-frame optical (based 
on the ground-based Near Infrared selected FIRES sam- 
ple) out to z ~ 2.5, and found that the average surface 
brightness at z ^ 2.5 is about 2-3 mag arcsec"^ brighter 
than in the local universe, but the average size at a fixed 
stellar mass has evolved by less than a factor of two. 



ITruiillo et al.l (|2006l . hereafter T06) presented the results 
of a similar analysis of a larger sample from FIRES, and 
combined those results with the lower redshift studies of 
S03 and B05. 

How do these observations compare with theoreti- 
cal expectations? In the simplest version of the Fall- 
Efstathiou picture, under the assumption that halo mass 
density profiles are singular isothermal spheres (SIS; 
p{r) cx 1/r^) and neglecting the self-gravity of the 
baryons, we expect the size of a galactic disk that forms 
within a dark matter halo to scale as rdisk oc Arvir, where 
A is the dimensionless spin parameter and rvir is the 
virial radius of the dark matter halo. N-body simulations 
have demonstrated that the spins of dark matter halos 
are not correlated with halo mass or most other proper- 
ties, and the distribu tion does not evolve with time (e.g. 
iBullock et al]|2001al) . If the stellar mass of the disk is 
a constant fraction of the halo virial mass, then, in this 
simple picture, the expectation is that the average size 
of galactic disks of a given stellar mass will evolve as rvir 
evolves for halos of a given virial mass. In the currently 
favored ACDM cosmology, this would imply a decrease 
of a factor of ~ 1.7 out to z = 1 and a factor of ~ 3 
out to z '--^ 3. This rdisk oc Arvir scaling (hereafter re- 
ferred to as the SIS model) has frequently bee n used in 
the hterature as a theoret ical baselin e fe.g. iMao et al.l 
fT998t iFerguson et al.ll20nl B05,T06). IMao et al.l (fl998h 
found that the SIS model was consistent with the size 
evolution of disks out to z ~ 1 compared with the data 
available at the time, but the samples were tiny, and the 
observational sele ction effects wer e not w ell characterized 
or accounted for. iFerguson et al.l (|2004( ) also found that 
the average rest-UV sizes of rest-UV selected galaxies at 
1-4 ^ 2 < 5 were consistent with the SIS model, but the 
connection of rest-UV luminosity with stellar mass (or 
halo mass) is quite uncertain. Most recently, B05 and 
T06 concluded that the predicted evolution in the SIS 
model is considerably stronger than the observed evolu- 
tion of the rest-optical sizes in their stellar-mass selected 
disk samples. 

However, the SIS model neglects several important fac- 
tors that are believed to play a role in determining the 
size of galactic disks forming in CDM halos. 1) The mass 
density profiles of CDM halos are not SIS, but have a uni- 
versal f orm (known as the N avarro-Frenk- White (NFW) 
profile; iNavarro et aLlfTQQTf ). characterized by the con- 
centration parameter Cvir- The concentration parame- 
ter quantifies the density of the halo on small {'^ kpc) 
scales relative to the virial radius, and has an impor- 
tant impact on the structural parameters of the resulting 
disk. There is a cor relation betwee n halo virial mass M^ir 
and concentration (Navarro et al. 1997), though with a 
significant scatter (jBuUock et al. ,2001b) . and this mean 
halo concentration-mass relation evolves with time, in 
the sense that halo s of a given niass we re less concen- 
trated in the past (|Bullock et al.l l2001bD . 2) The self- 
gravity of the baryonic material may modify the distri- 
bution of the dark matter as it becomes condensed in 
the central part of the halo ("adiabatic contraction") 3) 
Disks with low values of A and/or large baryonic-to-dark 
mass ratios may not have sufficient angular momentum 
to support a stable disk. These unstable disks may form 
a bar or a bulge, and might no longer be included in a 
sample of 'disk dominated' galaxies. 
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There are of course numerous other potential compli- 
cations in the process of the formation and evolution of 
galactic disks, which we do not consider here (though we 
discuss some of them in ^jS]) . Here we present the predic- 
tions of a model for disk formation that improves on the 
simple SIS model by incorporating NEW halo profiles, 
adiabatic contraction, and disk instability. This model 
is based on the f ormalism and basic ingredients pre - 
sentedj^jBlumcnthaLeTS] (|1986[ ). iFlores et all (|1993D . 
and lMo et all (11998 *. MMW98). Our main result is that 
the predictions of this improved model are compatible 
with the rather weak observed evolution of the disk stel- 
lar mass-size relation out to z 1 reported by BOS. We 
also extend these predictions out to higher redshift z ~ 3, 
and find acceptable agreement with the results reported 
by T06. 

We discuss the ingredients of our model in Sj^l give a 
brief summary of the observational data in ^ present 
our results in SjH and discuss our results and conclude in 
Jjll We assume the following values for the cosmological 
parameters: matter density Qm = 0.3, baryon density 
rib = 0.044, cosmological constant JIa = 0.70, Hubble 
parameter Hq = 70kms~^ Mpc~^, fluctuation amplitude 
(78 = 0.9, and a scale-free primordial power spectrum 
ris = 1. 

2. MODEL 

The fundamental hypothesis of our model, following 
the Fall-Efstathiou paradigm, is that galactic disks form 
within massive, extended dark matter halos. The struc- 
tural properties and multiplicity functions of dark mat- 
ter halos in a given cosmology can be robustly predicted 
from N-body simulations and analytic fitting formulae. 
Therefore, the main challenge is to relate the structural 
properties of dark matter halos to the observable proper- 
ties of disk galaxies. Making this connection is the goal 
of our model, which we now briefly present. Note that 
our approach, and the exposition here, a re closely base d 
on the formalism p resen ted in [Bl umcnth al et al.l (jl986f ). 
lElores et all ^93), and lMo eiTali ((1998, MMW98). 

We assume that a fraction fd = rud/M^i^ of the halo's 
virial mass is in the form of baryons that are able to 
cool and collapse, forming a disk with mass rud- We 
denote the angular momentum of the disk material as 
Jd, and assume that it is a fixed fraction of the halo 
angular momentum Jh and that specific angular mo- 
mentum is conserved when the baryons cool and col- 
lapse. We also assume that the disk is thin, in centrifugal 
balance, and has an exponential surface density profile 
I](r) = So exp (— r/r^). Here and Sq are the disk 
scalelength and central surface density, and are related 
to the disk mass through — 27rEo''d^- 

The angular momentum of the disk is: 



Jrf = 27r / Vc{r)J:{rydr. 



(1) 



If we were to assume that the initial dark matter den- 
sity profile is a singular isothermal sphere, and to neglect 
the self-gravity of the disk, then the rotation velocity 
Vc{r) is constant and equal to Vvir, and the angular mo- 
mentum of the disk reduces to: 



Using the definition of the spin parameter, A = 

A|£;,,|i/2G-iAfvir"^^^ (iPeebleslHtel) . where Eh is the 
total energy of the halo, we can write: 



rd 



(Jd/rrid) 

(A/Mvir) 



(3) 



We now define fj = (Jd/"^!^)/(A/A^vir) to be the ratio 
of the specific angular momenta of the disk and the halo. 

If we assume all particles to be on circular orbits, the 
total energy of a truncated singular isothermal sphere 
can be obtained from the virial theorem: 



Eh = — 



2^vir 2 

Inserting this expression in eqn. [31 we obtain: 
1 



(4) 



(5) 



If we assume that the spin parameter does not change 
with cosmic epoch, this leads to the expected scaling 
rd{z) cx Tviriz). We refer to this model as the 'SIS' model, 
reflecting the assumption that the initial halo profiles are 
singular isothermal spheres. 

A more realistic model would incorporate cosmologi- 
cally motivated Nayarro-F renk- White (NEW) halo pro- 
files (H avarro et al.l Il997f ) and should also include the 
effect of the self-gravity of the baryons on the predicted 
properties of the disk. If the gravitational effects of the 
disk were negligible, then the rotation curve would be the 
same as that of the unperturbed NEW halo, which rises 
to a maximum value Knax at a radius ~ 2rs , where Vs is 
the NEW scale radius = Tvir/cvir- However, on scales 
where the mass of baryonic material becomes significant 
compared with the dark matter, we must account not 
only for the direct gravitational effects of the baryonic 
disk material, but also for the fact that the dark mat- 
ter contracts in response to the gravitational force of the 
baryons as they collapse to form the disk. We compute 
the effect of this contraction under the assumption that 
the disk forms gradually enough that the response of the 
halo is "adiabatic", i.e., that the following "adiabatic in- 
variant" quantity is conserved: 



GMf{< rf)rf = GM,{< n)n 



(6) 



Jd = AnJ^oVcrd = 2rnd?'dKir 



(2) 



where Mi{< r^) is the initial mass within an initial ra- 
dius Ti (assumed to be given by the NEW profile), and 
Mf{< rf) is the post-collapse mass within a final ra- 
dius rj. Numerical tests of the validity of the "adiabatic 
collapse" formalism have shown that it works surpris- 
ingly well, even when some of the underlying assump- 
tions (such as spherical symmetry a nd all particles be- 
ing on circular orb its) are violated (jJesseit et al.ll200^ 
IGnedin et aLll20nl . 

The final mass within a radius r is the mass of the 
exponential disk plus the mass of dark matter within 
the initial radius r^: M/(r) = md{r) + M{ri){l — fd)- 
The modified rotation curve can then be written as a 
sum in quadrature of the contributions from the disk and 
the (contracted) dark matter halo: Vc^(r) = Vc^d^{r) + 
Vc,DA/^(f)- The rotation curve for the disk is given by 
the usual expression for a thin exponential disk (see e.g. 
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iBinnev fc TremaindllQSTl ). and the dark matter compo- 
nent of the rotation velocity is given by Vc^dm {r) = 
G[Mf{r) — md{r)\/r. Using eqn. [2] again, and making a 
change in variable u = r/r^;, we can write the expression 
for the angular momentum of the disk as: 

Jd = m,r,Kj e--u^^^4f^du, (7) 



(8) 



Jo ^ir 

and, in direct analogy to eqn. O we obtain 



rd 



where 



/fl(A,c,/d) = 2 



Vvir 



1 -1 



du 



(9) 



and fc = Etssfw /Esis is the total energy for a halo with 
an NFW profile, relative to that for a singular isothermal 
sphere with the same mass. Note that in this step we 
have also used the fact that rvir >> r^, and that the disk 
density profile declines exponentially, to set the upper 
limit of integration to infinity. 

For a given set of halo properties Mvir, Cvir, and A, 
and disk parameters fd and fj, we can solve this system 
of equations iteratively to obtain Vc{r) and rd We 
describe how we obtain the halo properties and disk pa- 
rameters in the following section. We will refer to the 
improved model as the 'NFW' model. 

Comparing eqn. [5] with the corresponding expression 
in the SIS model, eqn. [51 we can see that they differ by 

— 1/2 

two factors: fc , which reflects the difference in energy 
of a singular isothermal sphere vs. a NFW profile, and 
fn, which reflects both the NFW density profile and the 
effect of the adiabatic contraction. The ratio of eqn. [8] 
and eqn. [5] (rd/riso) is plotted in Fig.[Tl as a function of 
disk fraction fd and of halo concentration Cvir- From this 
plot, we can gain several insights into the behavior of the 
model. First, we can see that the NFW model generally 
predicts disk sizes that are smaller than the values in the 
SIS model. For a dark matter halo with Afvir = lO^^M©, 
the present-day virial radius is rvir — 260 kpc, and so 
from eqn El the predicted scale radius would be about 
6 kpc for the median value of the spin parameter A = 
0.035. This is almost a factor of two larger than one 
might expect, if such a halo hosts a galaxy similar to the 
Milky Way or M31. Therefore, there is room for this 
reduction in size predicted by the new model. 

We can also see that smaller values of fd lead to larger 
values of r^. This is because there is less mass and there- 
fore less gravity to contract the halo, leading to a more 
extended disk. Finally, we see that lower values of con- 
centration Cvir also lead to larger (more extended) disks. 
Lower values of Cvir imply that there is less mass near 
the center of the halo relative to the outer parts, again 
leading to weaker gravitational forces on the kpc scales 
where the disk is forming, and less contraction of the 
halo, therefore a more extended disk. 

For a more detailed account of how to go about solving 
the equations, and analytic fitting functions for fn and fc, sec 
MMW98) 



2.1. Disk Stability 

Because of the rather large range of allowed spin pa- 
rameters A, our model predicts a fairly broad distribution 
of galaxy sizes at fixed stellar mass. However, it is possi- 
ble that not all of these disks gravitationally stable. It is 
well-known from numerical studies that galaxies in which 
the self-gravity of the baryons domi nates may be unsta- 
ble t o the formation of a bar (e.g. IChristodoulou et al.l 
I1995D . Following MM W98, who base their st ability cri- 
terion on the study of lEfstathiou et all (|1982| ). we adopt 
a threshold for instability = T4iax/(Gmd/r£;)^/^ < 
em.crit, where V^nax is the maximum rotation velocity, 
which we approximate as the rotation velo city at ~ 3rd- 
We adopt em.crit = 0.75 (jSyer et al.lll999( ). and assume 
that galaxies with e„i less than this critical value may 
form a significant bar or bulge and be excluded from a 
disk sample. Because of the uncertainties associated with 
the fate of these "unstable" disks, however, we present 
the results of our analysis both including and excluding 
these objects. 

2.2. Initial Conditions and Model Parameters 

Let us consider a virialized dark matter halo of a given 
virial mass (Mvir) at a specific redshift. We require the 
following information, all as a function of redshift: 

• The number density of these halos 

• The halo virial radius and virial velocity 

• The initial halo mass density profile p(r), here char- 
acterized by the Navarro-Frenk- White concentra- 
tion parameter Cvir 

• The halo specific angular momentum, character- 
ized by the dimcnsionlcss spin parameter A. 

All of this information can be robustly obtained with 
a high degree of accuracy from modern numerical cos- 
mological N-body simulations. We now describe each 
ingredient in some detail, and how we obtain it. 
Halo number density as a function of redshift: Over a 
given redshift in terval, we select halo m asses from the 
mass function of lSheth fc TormenI ()1999l ) using a Monte 
Carlo procedure. 

Halo virial radius and virial velocity: It is standard 
to define dark matter halos as spherical regions within 
which the average overdensity exceeds a critical value 
Avir times the mean background density of the universe. 
The value of Avir is computed based on the spherical 
collapse model. In an Einstein-de Sitter (17 = 1) uni- 
verse, Avir = 178 and by definition the background den- 
sity pb equals the critical density Pcrit- During the early 
years of the development of the CDM paradigm, when 
the Einstein-de Sitter model was widely assumed to be 
correct, it became common to define halos as overdensi- 
ties greater than 200 times the critical density (the value 
200 was obtained by just rounding off 178). This as- 
sumption is still commonly used by many authors today, 
in spite of the fact that it is no longer well-motivated 
in the now-favored concordance cosmology. We instead 
use the appropriate value of Avir for our adopted cos- 
mology, as predicted by the spherical collapse model 
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Fig. 1. — The ratio of the disk scale length predicted in the NFW model with adiabatic contraction, to that in the SIS model 
(neglecting disk self-gravity), as a function of the fraction of baryons in the disk fd (left panel) and the NFW halo concentration 
Cvir- For all curves, the mean value of the spin parameter A = 0.05 has been assumed. In the left panel, the solid curves are for 
Cvir=10, while the upper and lower dashed curves are for Cvir=5 and 20, respectively. In the right panel, the solid curve is for 
fd ~ 0.06 and the upper and lower curves are for fd ~ 0.03 and fd = 0.09, respectively. 



(Avir(^ = 0) ~ 337) times the mean background den- 
sity. This assumption is also commonly used in the liter- 
ature, and was used in the Bullock et al. (2001a, 2001b) 
papers from which we obtain our halo profile and spin 
relations. The details of the adaptation of the spherical 
collapse model to a general cosmology (with cosmological 
constant) and formulae for and plots of the relations for 
halo virial radius as a fun ction of mass and redshift ar e 
given in the Appendix of ISomerville fc Primackl (|1999D . 
and are the same relations that we use here. The halo 
virial velocity is then easily obtained by use of the virial 
relation, Vviv = (GMvir/^'vir)^^^- Note that, because ha- 
los are defined with respect to the background density 
of the universe, average halo densities were higher in the 
past. This implies that halos of a given virial mass had 
smaller virial radii, and higher virial velocities, in the 
past than they do today. 

Halo Density Profile: We assume that initially the halo 
mass density profiles obey the Navarro- Frenk- White form 
(jNavarro et al.l [19971 ). and compute the halo concentra- 
tion Cvir as a function of mass and redshift usi ng the 
analytic model provided by B ullock et al.l ()2001bl . BOl), 
which is based on the results of numerical simulations. 
The results of BOl have been confirmed by a number of 
other studies (lJinell2000l: lEke et al.l l2001l: IWech sler et all 
[200llMacci6 et al.ll2006D . 

The basic ingredients of the BOl analytic model for 
Cvir(M, z) are as follows. They define a "collapse red- 
shift" Zc for each halo as M*(zc) = FM^ir, where 
fT[M*(z)] = Sc/D\in. Here, F is a parameter of the 
model, which is measured by fitting to N-body simu- 
lations, a{M) is the z = linear rms density fluctu- 
ation on the comoving scale encompassing a mass M, 
6c — 1.686 is the linear theory value of the critical den- 
sity in the spherical tophat model, and -Diin(z) is the 
linear growth factor. Then the concentration is given by 
Cvir(Afvir, 2) — K{1 + Zc)/(1 + z) , whcrc K is another 
adjustable parameter of the model. Note that Zc is in- 



dependent of z, and therefore Cvir oc (1 + z) ^ at fixed 
mass. 

iBullock et al] ()2001bD found values of = 0.01 and 
K = 4 for a ACDM model with the same values of 
Q\, and Hq as our cosmolo gy, but with a highe r value 
of (Tg = 1- A recent study bv lMaccio et al.l ()2006[ ). based 
on a larger sample of higher resolution simulations, con- 
firms the BOl model but finds K = 3.4 for a ACDM 
cosmology with as = 0.9 (corresponding to a 15 percent 
lower normalization than BOl, but with the same mass 
dependence). In the mod el presented here, we adopt 
the updated iMaccio et al.l (2006) normalization of the 
Cvir(-Mvir) relation {K — 3.4). 

The concentration of halos of a given mass at a given 
redshift depends on the cosmology and the power spec- 
trum. In cosmologies in which halos assemble their mass 
early (late), a halo of a given mass has a higher (lower) 
concentration riWechsler et al.l[2003 ). We can see from 
the above formulae that in a cosmology with a mod- 
ified D\in{z) (for example, if we changed the value of 
the cosmological constant, r^A), the redshift evolution 
of Cvir(-/Wvir) would bc different. Similarly, if we modi- 
fied the normalization or shape of the power spectrum 
(which enters through (t(M)), it would alter the normal- 
ization and slope of the Cvir(Afvir) relation. For example, 
in the cosmology derived from the WMAP three year re- 
sults of Spergel et al. (2006)^'^, the Cvir(Afvir) relation has 
a normalization about 30 % lower than in our adopted 
cosmology, but a nearly identical redshift dependence. 
Halo Specific Angular Momentum: N-body simulations 
have demonstrated that the spin parameter A, which 
characterizes the specific angular momentum of dark 
matter halos, is u ncorrelated with t he halo's mass and 
concentration (Bul lock et"ani2001at iMaccio et al]|2006l) 
and does not evolve with redshift. The distribution of A 

13 Q„ = 0.24, Qt = 0.042, Ha = 0.76, Ho = TSkms"! Mpc"!, 
(78 = 0.74, ria = 0.95 



6 



Somerville et al. 



is log-normal, with m ean A = 0.05 and width ax — 0.5 
(|Bullock et al.ll2001a^ . We therefore assign each halo a 
value of A by selecting values randomly from this distri- 
bution, assuming that it is not correlated with any other 
halo properties or with redshift. 

2.3. Predicting observable disk properties 

We assume that each halo with a virial velocity below 
350 km/s hosts one disk galaxy at its center. Once the 
halo properties are specified as described above (based on 
dissipationless N-body simulations), there are only two 
free parameters in our model: the fraction of mass in 
the form of baryons in the disk, fd, and the fraction of 
the specific angular momentum of the DM halo captured 
by the disk, fj. Throughout this work, we assume fj = 
1. Although there are physical reasons to think that fj 
may not be exactly equal to unity, note that in order 
to change our main conclusions, fj would have to be a 
strong function of redshift. There is no obvious physical 
reason to expect this to be the case. 

Since fd is now the only free parameter, we can fix its 
value by requiring that we reproduce the observed scale 
length of a galaxy of a given virial mass scale. For exam- 
ple, adopting fd — 0.06 produces a galaxy with disk mass 
rud = 6 X IG^^Mq and scale length = 3.6 kpc in a halo 
of mass Mvir = IO^^Mq, similar to the Milky Way and 
in good agreement with the median value for the scale 
length of disks in our SDSS sample at this stellar mass. 
However, if we assume a constant value of fd = 0.06 
for all halo masses, we find that the slope of the stellar 
mass vs. disk scale length relation is too steep, i.e. the 
sizes of smaller mass disks are too small, and those of 
larger mass disks too large, compared with the observed 
relation. This has also been noted be fore by other au- 
thors working with similar models (e.g. 'S hen et al]l2003t 
iDutton et aiT [2007). We can solve this problem by as- 
suming that the disk baryon fraction varies with halo 
mass, in the sense that lower mass halos form disks with 
lower baryon fractions. This is in fact what is found in 
physically motivated models of galaxy formation, as su- 
pernova feedbac k can m ore easil y heat and eject gas from 
low-mass halos (jNatarai an 1999: iDekel k Silklll986l ). We 
therefore adopt a simple functional form for fd as a func- 
tion of halo mass: fd = /o/[1.0 (Mvir/MJ"], following 
803, and find good agreement with the observed stellar 
mass vs. disk scale length relation with parameter values 
/o = 0.13, Afc = 1.0 X lO^^M©, and a = -0.67. Note 
that we do not discriminate between stellar mass and the 
mass in cold gas. For the relatively massive disks that 
we will focus on, the cold gas fraction should be fairly 
low. 

3. SUMMARY OF OBSERVATIONAL DATA 

For the main part of our analysis, we use the same 
sample of disk-dominated galaxies that was used for the 
analysis of B05. We give a brief summary of that sam- 
ple here, and refer to B05 for details. For our local 
z ^ sample, we use the NYU Va lue-added Galaxy Cat- 
alog (VAGC: [Blanton et al.ll200"5l ). based on the second 
data release of the SDSS (DR2). The VAGC catalog 
was used to obtain Sersic parameters, r-band half-light 
radii, and magnitudes. We then construct a sample of 
disk-dominated galaxies by requiring Sersic parameter 
n < 2.5. We convert the r-band half-light radii to rest- 



frame T^-band radii using the conversion derived in B05, 
Re{V) = 1.011i?e('')- We compute stellar masses from 
the (k-correcte d) q and r -band photometry, using the 
prescription of iBell et al.l (|2003f) . which relies on a con- 
version between g — r color and average stellar mass-to- 
light ratio. We a ssume a norma lization for this relation 
consistent with a lKroupal (|2001h IMF. 

Our high red shift dis k sample comes from the GEMS 
survey (jRix et a l. 2004; Caldwell et al]l2006f ). which con- 
sists of Vgoe a-nd zgso imaging over an area of 900 
arcmin^ with the ACS on HST. The 5cr point source 
detection limit is 28.3 magnitude in the Vgoe band and 
27.1 magnitude in zgso. The GEMS o bject catalog is 
based on the zgso image; for details see I Caldwell et al.l 
(|2006h . High -accuracy photometric redshift estimates 
((Tz/(z -f 1) ~ 0.02) are obtained from the ground-based 
COMBO-17 survey (iWolf et all 12001 . The i?-band se- 
lection limit of COMBO-17 {mn ^ 24) limits the range 
of the sample to redshifts z < 1. The 17-band pho- 
tometry of CO MBO has also been used to obtain stellar 
mass estimates ()Borch et al.ll2006[ ). assuming a Kroupa 
IMF. The GEMS main sample consists of almost 8000 
galaxies with COMBO counterparts and redshifts. We 
fit each galaxy with a Sersic profile and select a disk- 
dominated sample with good quality fits, Sersic n < 2.5, 
and extended light profiles. This sample contains 5664 
objects. Typical un certainties are 35% in r^ and ~ 0.2 
magnitudes in ()Haussler et al.|[2006[) . The apparent 
half-light sizes measured in the observed zs5o band are 
converted to rest-frame V-band using an average color 
gradient correction based on a local sample of disk galax- 
ies (see B05). These corrections are small (±3 %) over 
the entire redshift range of our sample (the observed zgso 
band samples the rest- frame y-band at z ~ 0.5). Where 
disk scale radii are quoted in this work, we have obtained 
them by simply assuming that the measured half-light ra- 
dius and the disk scale radius are related by the standard 
expression for a pure exponential disk, r^ = re/1.68. 

In order to estimate the completeness of the combined 
GEMS-I-COMBO disk sample, we have performed exten- 
sive simulations (Haussler et al. 2006; Rix ct al. 2004). 
Artificial disks were inserted into blank sky, and the 
source detection and fitting software was run on this 
image. Poor fits are excluded in the same manner 
as for the real galaxy images. We can then calcu- 
late the success rate for detecting and obtaining a good 
fit for the artificial galaxies, as a function of appar- 
ent effective radius and apparent magnitude. We mul- 
tiply this GEMS completeness factor by the (redshift, 
magnitude, and color-dependent) probability that the 
galaxy would be detected and successfully assigned a 
redshift in the COMBO survey. Based on these esti- 
mates, BOS argue that GEMS is not surface brightness 
limited even in the highest redshift bin, and that the 
combined GEMS-I-COMBO sample is complete down to 
stellar masses of lO^^M©. As in B05, we limit our anal- 
ysis to galaxies with stellar mass greater than this value, 
and we weigh galaxies by the inverse completeness factor 
in computing distributions and means. To avoid using 
galaxies with very large weights, we exclude objects with 
a completeness factor smaller than 0.5. B05 have shown 
that the average sizes and surface densities are insensi- 
tive to the choice of this limiting completeness factor at 
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our adopted stellar mass limit. 

4. RESULTS 

Using the models outlined above, we construct a mock 
catalog for a low redshift slice (0.001 < z < 0.2) repre- 
sentative of the SDSS sample, as well as a lightcone from 
0.1 < z < 1.1 with approximately ten times the area of 
GEMS (10 X 900 arcmin^), which we divide into five bins 
with roughly equal comoving volume (0.1 < z < 0.56, 
0.56 < z < 0.74, 0.74 < z < 0.87, 0.87 < z < 1.0, 
1.0 < z < 1.1). Following MMW98, we define "sta- 
ble disks" as those with stability para meter Cm. a.bove a 
critical value, and adopt em,crit = 0.75 (jSver et al.lll999D . 
The stellar mass vs. disk scale length relations for the six 
redshift bins from z ~ 0-1 are shown in Figure [21 for the 
mock catalogs on the left and the observed SDSS and 
GEMS disk samples on the right^'*. The observational 
samples are corrected for completeness as described in 
Section [3l 

4.1. The size-mass relation and its evolution to z ^ I 

In the left panels of Fig. [21 the contours show the re- 
sults for stable disks only. To illustrate the impact of ex- 
cluding the unstable disks, the median and 10 and 90th 
percentiles in size as a function of stellar mass for all 
disks, without any stability criterion applied, are shown 
in the z = 0.1 panel. The relative impact of applying 
the stability criterion on the size-mass relation and its 
scatter is similar in all the redshift bins considered here, 
so for clarity we show only the median relation for stable 
disks in the other panels. It is clear from the definition 
of the stability parameter that at a given halo mass, 
disks with larger values of fd are more likely to be unsta- 
ble, and also at a given stellar mass, more compact disks 
(low rd) are more likely to be unstable. Therefore, more 
massive galaxies are more likely to be unstable (recall 
that we assumed that fd increases with increasing halo 
mass), and so excluding the unstable (compact) disks re- 
sults in a steepening of the slope in the m^-r^ relation 
at m, ~ 2 — 3 X lO^^M©. Excluding the unstable disks 
also significantly reduces the scatter in disk size at fixed 
stellar mass for the massive disks. The fraction of disks 
deemed stable by our criterion is nearly constant over 
the redshift range considered here, and ranges from 95% 
at zr^ 0.1 to 97% at z - 1. 

Figure [H shows that the NFW model can reproduce 
both the observed slope of the size-mass relation and the 
scatter in size at a given stellar mass fairly well at z '--^ 0. 
Our results agree with those of 803 and other investiga- 
tions. The median and 10 and 90th percentiles in disk 
scale length as a function of stellar mass for the stable 
model disks from the low redshift bin are repeated in 
every panel, and against these contours we can immedi- 
ately see that the model predicts that the average size 
of disks at fixed stellar mass has increased by about 15- 
20% since z ~ 1. The GEMS sample is consistent with 
no evolution in size at fixed stellar mass. Note that we 

Note that, strictly speaking, the models predict the scale 
length of the baryonic mass in the disk, while the observed scale 
lengths are for the rest-frame V-band light. We do not attempt to 
correct for the known ~ 20% difference between these two quanti- 
ties, since we are mainly interested in the redshift evolution. We 
do note here, however, that time-evolving color gradients could 
therefore change our results. 



have normalized the model and data histograms to the 
same total number density in each redshift bin. 

An important side issue is that we could see in Fig- 
ure [2] that the model produces too many massive disks 
(m, > lO^^M©), especially at high redshift, compared 
with the observational samples. This is not surprising, 
since we have assumed that every dark matter halo con- 
tains a single disk galaxy, which is clearly not realistic, 
and we know that more massive halos (which produce the 
massive disks) are more likely to instead host an early 
type galaxy. Therefore, in order to avoid any bias from 
the unrealistically massive model galaxies, in Figures [3] 
andlH we show these distributions for galaxies with stel- 
lar masses in the range IO^^Mq < m* < IO^^Mq. 

Figure [3] shows the distribution of stellar surface den- 
sities for the same six redshift bins from z 0.1 
to z ^ 1. We define the stellar surface density as 
logS* = logm* — 2 log re — log(27r), where r^ = l.GSrd- 
From this figure, we can readily see that if we had not 
excluded unstable disks, the model would have predicted 
a broader distribution of surface densities than is seen 
in the data, with a highly skewed tail to very high den- 
sities that are not observed in the disk samples. With 
unstable disks excluded, the width of this distribution 
is reasonably consistent with the data at all redshifts, 
though there are hints of some interesting discrepancies; 
namely, the observational distributions appear perhaps 
a bit narrower and a bit more skewed than the model 
predictions. 

Figure [4] shows the evolution of the average value of 
logS, for the same stellar mass range. As seen in the 
previous figure, the model predicts significant, but fairly 
mild, evolution in log S, over the redshift range consid- 
ered (about 0.2 dex, or about a factor of 1.5; consistent 
with E» cx r~^). We also show the SIS scalings both 
for r2oo and for rvir (see Both of these predict 

much more dramatic redshift evolution than is observed: 
as found in B05, and as seen here, the results from the 
GEMS analysis are consistent with no evolution in the 
average value of over the redshift range 0.1 < z < 1.1. 

B05 found that the 2-tT error bars on the average 
values of log S obtained from bootstrap resampling are 
^ ±0.04-0.1 dex, but this certainly underestimates the 
possible systematic errors. For example, the stellar mass 
estimates may be systematically incorrect if galaxies have 
more bursty star formation histories at high redshift, or 
the size estimates could be systematically biased by the 
increasingly irregular morphologies of high redshift disks 
(our fitting simulations assume perfectly smooth galax- 
ies). We estimate the overall systematic uncertainty in 
size at fixed stellar mass to be ~ 30%, and show repre- 
sentative error bars refiecting this in Figure [31 The pre- 
diction of the NFW model is a significant improvement 
over the SIS scaling, and it is in acceptable agreement 
with the GEMS data to z ~ 1 within these estimated 
uncertainties. 

The surface density of disks is important because 
nearly all semi-analytic galaxy formation models assume 
that the efficiency of star formation is determined by 

Note that the evolution is somewhat more rapid when the 
definition rvir is used, because the virial overdensity used to define 
the halo, Avir, evolves with redshift while with the r2oo definition 
it remains constant. Because Avir is larger at higher redshift, the 
halos are smaller in radius. 
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Fig. 2. — The relationship between stellar mass and disk scale length for a "local" sample (0.001 < z < 0.2) and in five redshift 
bins of approximately equal comoving volume from z — 0.1 to z = 1.1. In the left panels, contours show the NFW model 
predictions for stable disks. In the right panels, contours show the completeness-corrected distributions for the SDSS and GEMS 
samples for the same redshift bins. The left and right panels are normalized to the same total number density in each redshift 
bin. The dark blue solid and dashed lines show the median and 10 and 90th percentiles, respectively, for the stable model disks 
in both columns of panels (models and data). The magenta solid and dashed lines in the left z = 0.1 panel show the median 
and 10 and 90th percentiles for all disks, without any stability criterion applied. The gray-blue dotted lines, repeated in each 
panel, show the z = 0.1 medians and 10 and 90 percentile lines for the stable model disks. 
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Fig. 3. — The distribution of stellar surface densities for disks with stellar mass 10'"Mq < m* < 10" M©, for a "local" 
sample (0.001 < z < 0.2), and in five redshift bins of approximately equal comoving volume from z = 0.1 to z — 1.1. Shaded 
histograms show the completeness-corrected distributions of SDSS and GEMS disks, selected via Sersic fits to their radial light 
profiles (n < 2.5). Solid squares show the mean value of log E* derived from these observations. Magenta dashed lines show the 
distributions for all model disks, and dark blue solid lines show the results for "stable" model disks only (see text). The solid 
dots indicate the mean of the model distribution for the current redshift bin, and the open dots show the means from all lower 
redshift bins. The upper set of dark blue dots are for stable disks, and the lower set of magenta dots are for all disks. 



10 



Somerville et al. 




2.0 



0.0 0.2 



0.4 0.6 O.f 

redshift 



.0 



o 
II 

N 



N 




Fig. 4. — The redshift evolution of the mean stellar surface 
density. Solid squares show the mean of log E, for SDSS and 
GEMS disks with 10^°Mq < m, < lO^M©. Open magenta 
circles show the results for all NFW model disks in this mass 
range, and solid blue dots show the results for stable model 
disks only. The dashed curve shows the evolution in E« that 
we would expect if disk size scaled like r2oo, and the solid 
curve shows the evolution for disk sizes that scale like rvir (as 
in the SIS model; see text). The SIS model predicts that disk 
surface densities were considerably higher at 2 ~ 1, in conflict 
with the data. 




Fig. 5. — The redshift evolution of the average size of disks 
with stellar masses greater than 3 x 10^'' Mq, relative to the 
average size of disks at z = 0.1. Square symbols with error 
bars show the observational estimates from TOG, obtained 
by combining the SDSS, GEMS, and FIRES datasets. Open 
(magenta) circles show the NFW model predictions for all 
disks. Solid (dark blue) dots show the NFW model predic- 
tions for stable disks only. The green dashed and solid curves 
show the scaling of r2oo and rvir (respectively) for dark mat- 
ter halos of fixed mass, as in the SIS model. The NFW model 
predicts more gradual size evolution, in better agreement with 
the observations than the SIS model. 



Fig. 6. — The redshift evolution of the circular velocity of 
disks with stellar mass 5-7 x IO'^^Mq, relative to the average 
circular velocity of such disks at z = 0.1. Open (magenta) 
circles show the model predictions for all disks. Solid (dark 
blue) dots show the model predictions for stable disks only. 
The green dashed and solid curves show the scaling of V200 
and Vvir (respectively) for dark matter halos of fixed mass, 
as in the SIS model. The NFW model predicts more gradual 
evolution in rotation velocity, in qualitative agreement with 
TuUy-Fisher observations to z ~ 1 (see text). 

the disk surface density through the Kennicutt Law, 
S* (X above a critical surface density Scrit- Pre- 

dicting disk surface densities correctly is therefore very 
important for self-consistent models of galaxy forma- 
tion. Note, however, that many semi-analytic rnodels 
(e.g.lSomerv ille fc Primacklll999l : lKauffmann eral]ll999l : 
ICroton eTa l. 2006) effectively assume that the disk sur- 
face density scales according to the SIS prediction. This 
will clearly lead to fairly serious inaccuracies in the pre- 
dictions of the evolution of the star formation rates in 
galaxies in these models. 

4.2. Evolution to higher redshift 

We now briefly explore how the models fare in com- 
parison with the more limited data available at higher 
redshift. Figure [5] shows the average size of disk galax- 
ies with stellar mass greater than 3 x lO^^Mo predicted 
by the model, out to z ~ 3, compared with the com- 
bined results from SDSS, GEMS, and FIRES, presented 
by T06. Each set of points is normalized relative to the 
average scale length of that sample at z = 0.1. Results 
are shown for all disks, and for "stable" disks only. Note 
that the absolute average disk scale lengths for the "sta- 
ble" disks are always larger than those for the total set 
of disks, with no stability criterion applied, however, the 
redshift evolution of the stable disk sample is slightly 
steeper than that of the overall sample out to 2; ~ 2. 
This is why the stable disk sizes, when normalized rel- 
ative to the average size of stable disks at z = 0, are a 
little bit lower than the total disk sample out to z ~ 2. 
We again also show the scaling for halo virial radius at 
fixed halo mass, which would predict much more rapid 
size evolution than is observed. As we have already seen, 
the improved model predicts fairly mild evolution in the 
average disk sizes out to z '--^ 1 , in quite good agreement 
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with the GEMS data. At higher redshifts z ~ 2-3, our 
model predicts that disks should be about 60% as large 
as they are today at a given stellar mass. This represents 
somewhat more evolution than the observational results 
of T06 indicate, but is within the quoted error bars. 

4.3. Evolution of the Tully-Fisher relation 

Our models also simultaneously predict the rotation 
curves of the disks. It has been shown in previous work 
that models similar to the one presented here can re- 
produce the d etailed rotation curves for the Milky Way 
and M31 (Kl vpin et al.ll2002f ). and the Tully-Fisher re- 
lation (circular velocity at some fiducial radius vs. lu- 
mino sity, or vs . ste l lar or baryon i c ma. s s) at redshift 
zero (|Mo et al.l 119981 : iDutton et all l2007t iG^edin et all 
[20061) . A detailed comparison with observations of ro- 
tation velocities and the Tully-Fisher relation is beyond 
the scope of this paper, but we briefly present a pre- 
diction for the evolution of the Tully-Fisher zero-point 
(here expressed as the average rotation velocity at ^3 
disk scale lengths, for disk galaxies of a fixed stellar 
mass, m* ~ 6 X IO^^Mq) in Fig. El We find that the 
disk rotation velocity Vdisk = Vc{r = Sr^) at fixed disk 
mass evolves much more gradually than the halo virial 
velocity Kir at fixed halo mass. In fact, our model pre- 
dicts that Vdisk should actually remain nearly constant 
from z ^ 1.5-3. Observational constraints on the time 
evolution of the Tully-Fisher relation, particularly in 
terms of stellar or bar yonic mass, remain uncertain, but 
IConselice et al.l (|2005l l found no evolution i n the K-band 
or ste llar mass Tully-Fisher to z ^ 1, and iKassin et al.l 
()2007t l also find no evolution to z ^ 1 in a redefined 
version of the stellar mass Tully-Fisher relation^^. 

5. DISCUSSION AND CONCLUSIONS 

We have shown that a CDM-based model of disk for- 
mation produces good agreement with the observed weak 
redshift evolution of the disk size-stellar mass relation 
from GEMS out to z ~ 1. This result is in contrast to the 
considerably stronger evolution implied by the assump- 
tion often used in the literature, that disk sizes simply 
scale in proportion to the dark matter halo virial radii. 
Similarly, we find much weaker evolution in the disk rota- 
tion velocity at a given disk mass than we would expect 
if Vdisk scaled in proportion to the halo virial velocity 

Vvir . 

The more gradual evolution in the "improved" model 
results primarily from the use of realistic NFW halo pro- 
files, and from the incorporation of the redshift evolution 
of the halo concentration-mass relation predicted by N- 
body simulations. Recall that in our model, at fixed halo 
mass, there are two halo properties and two disk parame- 
ters that determine the structural properties of the disk: 
halo concentration Cvir, halo spin A, the baryon fraction 
of the disk fd and the fraction of specific angular momen- 
tum captured by the disk materal, fj. We have assumed 
that the distribution of spin parameters A for the overall 
population of dark matter halos does not change with 
time, as demonstrated by N-body simulations. We have 
also assumed that fd is a fixed function of halo mass 

Instead of using the usual rotation velocity, they define a 
new velocity variable which includes contributions from random 
motions as well as rotation. 



that does not depend on redshift (though in reality this 
may not be true), and we have assumed a fixed value of 
fj — 1. Therefore, there is only one ingredient left that 
can possibly be responsible for changing the predictions 
for the redshift dependence of disk sizes at fixed mass: 
the halo concentration vs. mass relation. 

As we have discussed, simulations have shown that the 
halo concentration vs. mass relation depends on redshift: 
the h alo concentration at fi xed mass scales as Cvir oc (1 -I- 
(|Bullock et al.ll2001b[ ). Therefore, a halo of a given 
mass is less concentrated at high redshift. This apparent 
evolution, however, is really a consequence of the way 
that halos are assembled in a CDM universe. Studies 
of the mass accretion history of halos in simulations has 
shown that they have two basic phases of growth: an 
early, rapid phase, in which the central de nsity is set, and 
a seco nd phase of more gradual accretion (jWechsler et al.l 
I2OOI) . The mass within the characteristic scale radius 
rs is assembled during the early, rapid accretion phase. 
Afterwards, Ts stays nearly constant, while rvir increases 
due to smooth accretion of mass, leading to a formal 
decrease in Cvir = rvh-ffs- 

In Section [2] we demonstrated that less concentrated 
halos produce larger disks because of the lower mass and 
weaker gravitational forces in the central parts of the 
halo, where the disk forms. Thus the collapse of baryons 
produces less contraction of the dark matter profile, and 
a more extended disk. The trend towards lower con- 
centrations at earlier times therefore counteracts the de- 
creasing virial radii. Out to about z ~ 1, these compet- 
ing effects nearly cancel out, leading to the weak evolu- 
tion in the size-mass relation that we have shown. 

Note that we have repeated our calculations us- 
ing the "WMAP3" cosmological parameters derived by 
ISpergel et al.l (|2006D . Because of the reduced small scale 
power in this cosmology, halos form somewhat later and 
therefore have lower concentrations for their mass at the 
present day. However, once we renormalize our model 
(by making a minor adjustment to the normalization 
of the disk fraction relation fd) to again reproduce the 
size-mass relation of present day disks, we find that the 
redshift evolution of disk sizes and circular velocities is 
nearly unchanged. This is not surprising, as the redshift 
evolution of the halo concentration vs. mass relation is 
nearly identical to the one we originally adopted. 

We see a hint that the evolution predicted by these 
models is still a bit stronger than that indicated by the 
data. This could be a sign that one of the other assump- 
tions in our simple model is incorrect. For example, if 
the disk baryon fraction /d at a given halo mass decreases 
with increasing redshift, this would lead to shallower evo- 
lution and relatively larger disks at high redshift. Be- 
cause we have measured the disk sizes in the rest-V^ band, 
evolving color gradients could also mask evolution in the 
true size of the stellar disk. Alternatively, systematic 
biases in our stellar mass and size estimates could be 
impacting the observational estimates. 

We also see an increasing level of discrepancy at higher 
redshifts, z > 1.5. This could be a hint that an entirely 
different mechanism could be responsible for setting the 
sizes of disks at very high redshift. Mergers between gas- 
rich disks could result in a iiew, more spatially extended - 
disk (iKazantzidis et aI.ll2005HSpringel fc Hernailis^l2005l : 
[Robertson et al.ll2006f) . This same scenario could also 
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help to explain the kinematics of Damped Lyman-a sys- 
tems, which are difficult to reconcile with the standard 
Fall-E f st at hiou picture of disk formation (|Maller et al.l 
[2OOI . 

We also presented a prediction for the redshift evolu- 
tion of the average rotation velocity at ~ 3rd (close to the 
maximum of the rotation curve) for disks of fixed stellar 
mass TO* ^ 6 X 10^°Mq (i.e., the zeropoint of the TuUy- 
Fisher relation). We found a similar result as we did 
for the disk sizes: Vdisk decreases with time, but much 
more slowly than we would expect if we assumed that 
Vdisk « Vvir- The physical reason for this is the same as 
the one we have discussed in relation to the weak size 
evolution: halos have lower concentrations at high red- 
shift, and thus less mass near the center and so lower 
rotation speeds at the radii where rotation velocities are 
measured (near the maximum value of the rotation curve, 
Knax)-Thus Fmax even for pure DM NEW halos, with- 
out accounting for baryo ns, evolves more slowly than 
Vvir ()Bullock et al.|[2001bf ). In addition, in our model, 
lower concentrations lead to less contraction and more 
extended disks. This leads to flatter, less "peaky" rota- 
tion curves, and lo wer values of V d isk relative to a mo re 
concentrated halo ()Mo et aLlUQQllKlvpin et al.ll2002f) . 

About the same fraction of disks are classified as un- 
stable according to the condition we adopted (e™ < 0.75) 
over the whole redshift interval < z < 3, and the ex- 
clusion of unstable disks from the sample changes the 
average size by nearly the same amount over this inter- 
val as well. Therefore, as implemented here, disk stability 
does not play a significant role in determining the relative 
time evolution of the stellar mass-size relation, although 
it is important for reproducing the correct distribution 
of disk stellar surface densities. At very high redshifts, 
z > 2.5, our model predicts that the number of unstable 
disks starts to increase. However, we do not even know 
that thin, rotationally supported disks exist at such high 
redshifts, so we do not know how seriously to take this 
prediction. 

While the model presented here represents a significant 
improvement with respect to the overly simplified Arvir 
scaling commonly used in the literature, it still neglects 
many important aspects of disk formation in a hierar- 
chical universe, in particular the impact of mergers. We 
have also ignored the possible presence of spheroids and 
cold gas in our disk galaxies. As well, the fraction of 
baryons in the disk component as a function of halo mass 
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